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f~>) ' Abstract. We study the Glauber dynamics for the zero-temperature stochastic Ising 

I model in dimension d ^ 4 with "plus" boundary condition. Let 7+ be the time needed 

for an hypercube of size L entirely filled with "minus" spins to become entirely "plus" . 
O ' We prove that 7+ is 0{L^ {log L^) for some constant c, not depending on the dimension. 

, This brings further rigorous justification for the so-called "Lifshitz law" 7+ = 0{L^) 

[ni E] conjectured on heuristic grounds. The key point of our proof is to use the detailed 
knowledge that we have on the three-dimensional problem: results for fluctuation of 
monotone interfaces at equilibrium and mixing time for monotone interfaces dynamics 
(~| ' extracted from [2], to get the result in higher dimension. 
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1. Introduction 

^ ' We study the zero-temperature stochastic Ising model in a finite hypercube of side- 

\Q . length L in with plus boundary condition. Initially, the whole cube is filled with minus 

\ spins. The spins evolve following the majority rule: with rate one each spin takes the same 

value as the majority of its neighbors (including those at the boundary) if the latter is 
well defined. Otherwise its sign is determined by tossing a fair coin. Eventually, all minus 
spins disappear, as a result of the pressure imposed by plus spins from the boundary. 

Our aim is to study the time 7+ needed for the spins in the hypercube to become 
entirely plus. According to a heuristics from Lifshitz on a macroscopic scale, each point 
of the rescaled interface between plus and minus spins should move feeling a local drift 
proportional to its local mean curvature (this conjecture has been recently proved with co- 
authors for the zero-temperature two-dimensional case with macroscopically convex initial 
domain of — spins: see |5J). This readily implies that, with high probability, 7+ = 0{L^). 

Here we prove that, with high probability, 7+ = 0(L^(log L)^'') for all dimension d ^ 4. 
Our result complements existing analogous bounds for dimension d = 2 and d = 3 obtained 
in [U [2]. It is important to keep in mind that the cases d = 2 and d = 3 are completely 
different, mainly because the equilibrium fluctuations of Ising interfaces in d = 2 and in 
d = 3 occur on very different scales {0{L^/^) for d = 2 and 0(log(L)) for d = 3). As 
such, they have been analyzed by using very different approaches. The case d ^ 4 should 
be more similar to the case d = 3, which therefore plays the role of a critical dimension. 
However, contrary to what has been done in [2], our proof does not attempt to control 
the local mean drift of the interface. Such an approach, in fact, would at least require two 
main missing tools: (i) a detailed analysis of the local equilibrium fluctuations of Ising 
interface^; (ii) good estimates on mixing times for (d — l)-monotone surfaces, following 
e.g. the method developed by Wilson [9j for the three-dimensional case. Instead, we use 
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-'^In dimension d J5 4 the fluctuations of Ising hyper-surfaces are believed to be order one. Recently a 
result of this sort has been proved under a Lifshitz condition [S]. 
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an ad-hoc construction that allows us to bound the hitting time 7+ for the d-dimensional 
dynamics using known sharp estimates for the same hitting time in three dimensions 
[2]. One of the reasons why this approach is successful is that in dimension three, Ising 
interfaces are flat with only logarithmic fluctuations (see e.g. Proposition 4 in [2] for more 
details) . 

Lifshitz's law is a general prediction for the mixing-time of the stochastic Ising model in 
the whole low-temperature phase and not only for the zero-temperature version presented 
here. However, rigorous results available at positive temperature (especially concerning 
the upper-bound) are still very far from it. We refer to [7j for the best available result 
in dimension 2 and to its bibliography for a full account on what is known at positive 
temperature. Getting polynomial upper-bound (in L)on the mixing-time in dimension 2 
and even sub-exponential bounds in dimension 3 are challenging open problems. Note also 
that our result can be extended to some cases where the temperature tends to zero as the 
size of the system L goes to infinity, as in [2]. 

2. Model and result 

2.1. Definition of the model. Given a finite subset T C Z'^ we study a continuous time 
Markov chain on the state space 

fir :={-!, +1}^. (2.1) 

This Markov chain can be seen as heat-bath Glauber dynamics for the Ising model at zero 
temperature. We write a = {crx)xer for a generic element of fir- We refer to ax = ±1 as 
spins, and will often write simply + or — . Let 

ar:={y GZ'^\r I 3xEr,x~y} (2.2) 

denote the boundary (in Z*^) of F and fix some r/ G { + 1, —1}'^'" (referred to as boundary 
condition). We define also the internal boundary of a set to be 

d-T := {y eV \3x eZ'^\r,x y}. (2.3) 

Given ^ G fir, we consider the dynamics (f''^(0)t>o starting from initial condition 
£7^(0) = ^ (one may also write (T^'^(i) when one wants to underline dependence with 
respect to the boundary condition) and with the following evolution rules: sites x G F are 
equipped with independent Poisson clocks with rate 1, and the spin at x may change its 
value only when the clock at x rings. If the clock at x rings at time to, then one looks at 
the spins of the 2d neighbors of x (some of which may lie in dT, in that case their spin is 
given by r]) just before Iq, and then 

• if a strict majority of neighboring spins are + then set (7x{to) = +, 

• if a strict majority of neighboring spins are — then set o'x{tQ) = —, 

• if there are exactly d pluses and d minuses in the neighborhood of x then set 
crtf^{x) = + or — with probability 1/2 independently of the rest of the process. 

Remark 2.1. One gets a natural construction of the dynamics by equipping every x G 
F with independent Poisson clock processes {Tn,x)n ^ o of rate 1 and with independent 
sequences {Xn,x)n^ a of i.i.d. Bernoulli variables of parameter 1/2 (taking values in 
{ — 1,+1}). Then the evolution is set as follows: 

At time r„,x5 if a strict majority of the neighbors of x have spin + (resp. — ) one sets 

= +5 (resp. — ), otherwise one sets c7x{Tn,x) = ^n,a;- 
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This construction (sometimes referred to as the graphical construction in the Uterature) 
is used in the proof, and it is very useful to provide couphng for stochastic comparisons. 

More formally, one can write explicitly the transition rates p for our Markov chain as 
follows. Given a G Qr and x G T let (T^'^,(resp. a^'~) denote the configuration which is 
identical to a on $7 \ {x} and where ax'^ = + (resp. ax'~ = —)■ Then 

1 if J^yr^x, yer + Ej/~x, year ""ly < 0' 

p(f7^'+, CT^'-) = \ 1/2 if Y.y^x, s/GF ^^y + Ey^x, yedF % = 0, (2.4) 

if J2yr^x, yer ^y + Ej/~x, year % > 0- 




if Ylyr^x, yer ^y + Ej/^x, year % < 0, 
p(^-.-,a-.+) = <j 1/2 if E,~.,,er^s/ + Ej;~.,,GOr% = 0, (2.5) 

if Ey~x, s/gf "^s/ + Ej/~x, yedv % > 0- 
All other transitions have rate zero. We denote by P the probability associated with this 
Markov chain (although boundary condition and domain considered may vary) and by 
E the corresponding expectation. We are interested more specifically in the case where 
r = := {1, . . . is the d— dimensional hypercube, and ry is uniformly + on OTl. 
In that case, the dynamics possesses a unique absorbing state: the configuration with + 
spins everywhere otiVi (for simplicity we write it +). We are interested in estimating the 
hitting time of this absorbing state starting from the all — configuration 

r+ := inf{t ^ I CT-(t) = +}. (2.6) 

In fact our main result is an upper bound on the mixing time defined here as 

T^ix := inf{t ^ | P[r+ > t] ^ 1/4}. (2.7) 

It has been conjectured that the growth of Tmix should be given by Lifshitz's law: ^Z^ix 
should be of order (see [6] and [3] where a more general formulation of this conjecture 
can be found). This result has already been proved in dimension 2 (exactly) and 3 (with 
non-matching logarithmic corrections for the lower and upper bound) by Pontes et al. [l] 
(exact upper bound in dimension 2) by Caputo et al. \2\ (lower bound in dimension 2 and 
result in dimension 3). In [5], the first order asymptotic of T^\x (Ilnix = I? jUX + o(l))) 
is given for d = 2. We refer to [2] and its bibliography for a more complete introduction 
to this problem and underlying issues. The aim of this note is to extend this result by 
showing that there exists some constant cq > independent of the dimension d such that 

T^^ = 0{L\\ogLf^). (2.8) 

Previously in [31 Theorem 1.3], the bound T^\^ = 0{L'^) had been derived from an iterative 
procedure using as a starting point the result for d = 2. The same kind of procedure 
starting from the result of [2j would give Tmix — 0(^L^ "'^ (log Z^)*^) which can be considered 
as the best bound that had been known so far (see Proposition 13.11 for the statement and 
its proof for d = 4) . 

Our proof strongly relies on estimates that have been obtained for the three-dimensional 
model [2]. Sadly enough, we cannot yet give a lower bound on Tmix with the same poly- 
nomial order. The reason for this is that methods relying on exact computation for d = 2 
and 3 developed in a paper of Wilson [9] and used in [2] fail to work in dimension larger 
4. 



We state now our main result in full detail. 
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Theorem 2.2. Consider the dynamics in '■= {1, • • • ,L}'^, with boundary condition + 
on dTi- There exists cq such that for all d ^ 4, the following holds for every k > 0. 
Uniformly for all initial configuration ^, 



BxGTl, 4{L\logL) 



(2.9) 



We stress once again that the constant cq in the result does not depend on the dimension. 
In fact, examining the proof in [2] a (non-optimal) choice can be to choose any cq > 19/2. 
In what follows, we fix C2 := 3/2, ci > 13/2 and cq := ci + 2c2 + £ (for some small e). 



Remark 2.3 (about the lower-bound). The best lower-bound available for 7+ in dimen- 
sion d ^ 4 is rather trivial: there exists a constant c^ such that, starting from the all — 
configuration, 

lim P[r+^ CdL] = 1. (2.10) 

To obtain such a bound, it is sufficient to consider that when 7+ occurs, the center of the 
cube has turned from — to + so that there has been a sequence of consecutive updates 
on a paths of nearest neighbor sites, starting from the boundary of the cube and ending 
in its center. Any improvement on this bound would be of interest. 



3. Preparatory work 

3.1. Monotonicity of the dynamics. We introduce a partial order on Or (resp. ^or), 
saying that a ^ a' ii ^ cr'^ for every x £ T (resp. x G dT). Our dynamics enjoys nice 
monotonicity properties with respect to this order. Let /if ^ denote the law of cr^{t) with 
boundary condition r]. Then if ry ^ ry' and £, ^ 

where ^ denotes stochastic domination (i.e. fi ^ v ii for any increasing function / : Or — ^ 
R, fi{f) ^ i^{f), a function / being increasing if f{a) ^ /(o"') for all a ^ a'). This can be 
proved by constructing a coupling of the two processes cT^'^(t) and (t), building the 
two processes as explained in Remark 12. H with the same realization of Tn^x and Xn^x for 
both. With this coupling one has that for every t ^ 

a^'^^it) ^ a«'''''(t). (3.2) 



This result remains valid with boundary conditions that vary through time: if r/(-) and 
r/'(-) are (deterministic or random) functions from M+ to Or such that r]{s) ^ r]'{s) for all 
s G [0,t], one can still construct the Markov chains cr^'^^'\t), a^'^^'\t) (in this case the 
transition rates are time dependent), and couple them in a way that (13. 2p holds. 

A simple consequence of this monotonicity property is that to prove our main theorem, 
we just need to prove it with initial configuration all — . This monotonicity of the dynamics 
is crucial for our proof as it allows us to use an explicit strategy to control the evolution 
of the set of — (see Proposition 13.21 below). 
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3.2. Main line of the proof. Prom now on, we focus on the case d = 4. We show at the 
end of the paper how to adapt the proof for arbitrary dimension d. 

The first tool of the proof is to perform surgical changes on the original dynamics that 
cancels some of the spin- flips from — to +. Due to the above mentioned monotonicity 
properties, this slows down the dynamics and makes the mixing time longer (which is 
alright to get an upper-bound) , but if these changes are done carefully, this can also guide 
the dynamic in a nice pattern which allows us to have a better control on it. 

We choose carefully our modifications so that the dynamics we obtained is mostly a 
product of three-dimensional dynamics. We consider our four-dimensional cube of side- 
length L (we will rather consider a cylinder in what follows), as a superposition of L 
three dimensional slices, and use our update blocking strategy to make the evolution of 
these slices independent. Then we use some ingredients borrowed from [2] to control the 
evolution of each slice. In order to do this without slowing down the dynamics too much 
one needs an ad-hoc geometrical construction that we present at the end of this section. In 
order to make the reader more aware of this need, one exposes first a naive and non-optimal 
strategy that yields consequently a non-optimal bound. 

3.3. A non-optimal strategy. A simple way to control the dynamics in to make the 
successive three-dimensional layers of our four-dimensional cube evolve one after the other. 
This gives the following result. 

Proposition 3.1. When d = A, there exists a constant c such that for every k ^ 0, 

P[n>L^ (log Lr]=0{L-'^) (3.3) 

Proof. One consider a a modified version of the dynamics a~ (starting from all — ), that 
blocks the update of spins in the i-th layer {1, . . . , L}^ x {i} up to time {i — l)L^(log L)'^. 
By monotonicity a is dominated by a and thus it is sufficient to prove the proposition 
with a replaced by a. 

We prove by induction on i = 1, . . . ,L that for L large enough, 

P[3xG {!,..., L}3x{l,...,i},5,(iL2(iogL)-) = _] <^ iL-''-\ (3.4) 

For 1 = 1, notice that only the sites in the first layer, {1, . . . ,L}^ x {1} can be updated 
during the time window [0, -L^(log L)'^]. For each x in this layer, spins oi x + and 
X — 64 (where 64 := (0,0,0,1)) are respectively — and -|- so that neighbors along the 
fourth dimension have no infiuence on the majority rule. As a result the dynamics in 
{1, . . . ,L}^ X {1} is in bijection with the three-dimensional dynamics on a cube with -|- 
boundary condition and by ^ Theorem 1] (strictly speaking the result says that the 
probability that a — spin remains in the cube after time L^(logL)'^ is 0{L~^) but from 
the proof it can be shown that it is 0{L~'^) for all k), (|3.4p holds with d= 1. 

For the induction step, notice that if all spins of {1, ... , L}^ x {1, . . . , i} are equal to -|- 
at a given time, they stay -|- forever (they will alway keep a strict majority -|- neighbors). 

Thus if a{iLp'{[ogLY) = + for all x G {1, . . . ,L}^ x {1, . . . ,i}, then during the time 
window [iLp'ilogLY, [i + l)L^(log L)'^], only the spins in the layer {1, . . . ,L}^ x {i + 1}, 
may change sign (those above are blocked by definition). Thus for the same reason as 
above, the dynamic in this layer is in bijection with the three-dimensional dynamics, and 
thus at time [i + l)L'^(log L)*^ if L is large (using the Markov property for the dynamics 
and [21 Theorem 1]) all the spins have flipped to -|- with probability at least L~^~^. 
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P [3x G {1, . . . , L}3 X {1, . . . , i + 1}, a^{{i + l)L2(log Lf) = -] 

^P[3xe{l,...,L}3x{l,...,i}, aSL\^ogLY) = -] + L-^-\ (3.5) 
which combined with the induction hypothesis, ends the proof. □ 

3.4. The geometrical construction. The main reason why the strategy of the previous 
section does not give a sharp bound is that the evolution of the modified dynamics is very 
far from the predicted motion by mean curvature. The whole challenge is thus to make 
the dynamics evolves in a way that is not too far from motion by mean curvature, but 
that we can still control. This is the object of the following construction: 

For any r ^ let Sf denote the three-dimensional discrete ball with radius r (r need 
not be an integer) 

53 := {z = {zi,Z2, Z3) G Z3 I + + z| ^ r^} . (3.6) 
We consider dynamics in the cylinder 

■= 5'|i(iogL)C2 X {1, . . . ,L}, (3.7) 

with + boundary condition at the bottom of the cylinder and around it and — on the top 
of it, or formally: 

• f?x = + for X G ^'S'|j;^(iog L)C2 ^ {!)••• L}, and 5£ x {0} where d denotes here the 
boundary in Z^, 

• Vx = - for 5|^(i„gi)c2 X {L + 1}. 

We call this boundary condition rjQ. We put — spins at the top of the cylinder only because 
it turns out to be handy in the proof but this is not crucial point. 

Note that from this choice of boundary condition and from the way we prove of Theorem 
12.21 (see below), that the conclusion of Theorem l2.2l still holds also if the boundary condition 
is not all + on the boundary of the cube, but is — one one side and + on all the others. 

In most cases, we do not underline dependence in the boundary condition and write 
{iji{t))xeCL for the spin configuration at time t starting from configuration ^ with rjQ 
boundary condition. 

For every i. A; G N one defines the sets SL,k and as follows : 

5L,fc := 48L-fe)(logL)-2, ^^^^ 
Cf> := Cl n {z = {ZI,Z2,Z3,Z4,) G I {zi,Z2,Zs) G 5i,(i_2^4+2)+} ) 

where x+ = max(2;,0) denotes the positive part of x. Note that C^^^ = Cl- The set C^*^ is 
a pile (along the fourth dimension) of L three-dimensional balls of increasing radius (see 
Figure [T]) . For different subsets T> G Cl we may consider 

vn{i:'x{j}) (3.9) 

the j-th slice of T>. Our main result is a consequence of the following proposition, that 
controls the evolution of the set of minus in cr~{t). We prove it by considering an auxiliary 

dynamics that roughly blocks the spins in C'l^ to - up to time {i - l)L(log L)"'' (in fact 
one needs to perform several coupling but this is the main idea). 
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Proposition 3.2. For every k > 0, for every i £ {1, . . . , 4L} one has 

Phte [iL{\ogLy'\L^], 3xECl\4'\ a-{t) = -] =0{L-^). (3.10) 

Remark 3.3. In addition to an upper bound for the mixing time, the above result gives 
a control on how the set of — spins reduces through time: with high probability, at time 
iL{\og LY" and after, the set of minus spins is included in C^^ . However, the real evolution 
of the shape of the set of — should be very different from the one described in Figure [TJ 
Recall that at a macroscopic level, it is believed that the local drift of the interface between 
+ and — is proportional to local mean curvature (this conjecture is the heuristic support 
for Lifshitz's law). 

We show now how to get our main result from the former statement, which we prove in 
the next section. 

Proof of Theorem \2. S\ from Proposition \3.2\ when d = 4. Applying Proposition 13.21 for i = 
4L, we get that for any given k, for L large enough and for t = tL := 4L^(logL)'^o 



P[3xeCL\{SlL^i,^Ly,x{l,...,L}), a-itL) = -\ =0{L-'). (3.11) 

In particular if L is chosen large enough, at time ti there are no more minus spins in 
the cube 

Vl := { r6L(log LY^ + 1, . . . , r6L(log Lf^} + L} x {1, . . . , if, (3.12) 

with probability at least 1 — 0{L~^) since C {Cl \ C^l^^) (Xl is a translation of 
from Theorem I2.2p . Moreover, by monotonicity in the boundary condition, the evolution 
of the spins vnTi for the dynamics in starting for all — is stochastically dominated 
by the evolution oiVi started from all — inside with + boundary condition on dV^, and 
therefore holds. □ 



3.5. The special role of dimension three. The proof of Proposition l3.2l that we develop 
in the next section, relies heavily on the work done for the three-dimensional model. Before 
starting it, we wish to discuss the particular role played by the three-dimensional model. 

In [2], the bound Tmix = 0(L^(log L)^") -^^s obtained by a careful analysis of the 
fluctuations at equilibrium and the mixing time for monotone surfaces. A monotone 
surface is the graph of a function taking integer values, defined on a connected set A dl? 
and that is monotone with respect to both coordinate in 1? (this notion can be extended 
to any dimension). 

The reason why dimension 3 allows to get the result for higher dimension is that fluctu- 
ation of monotone surface have low amplitude (order logL for surfaces of side length L), 
and the reason is the following: if we want the dynamic on the different slices of not 
to interact between one another, one want to put the +/— interfaces of successive slices 
at a distance strictly larger than the typical interface fluctuation. This is the reason for 
the particular construction of the sets . This forces us to consider dynamics on a set 
Cl that has diameter L(logL)'^2 instead of L, but the consequence of this is just the loss 
of some extra logL factors in the result. 

In dimension 2, fluctuations of monotone interfaces are too large (of order y/L for 
surfaces of side length L) to make a similar construction efficient (in particular one could 
not have proved the main result in [2] for d = 3 by using only estimates for the two 
dimensional model). In dimension 4 and more, fluctuations at equilibrium of monotone 
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8L(log ly^ 



c 




c 



(2L) 




6L(logL)^2 




4L(log LY^ 



Figure 1. Schematic two-dimensional representation of the sets for four values of 
i. The outside rectangle in full line is Cl, the region within the dotted line is C^'"*. The 
figure also shows how we place the hypercube Tl (shadowed region) in C[''^'\Ci for the 
proof of the main Theorem. For graphical reason the scaling in the different directions 
is not the same. 



hyper-surfaces are believed to be smaller (i.e. analogy with Gaussian Free Field, one 
predicts that the mean fluctuations are of order 1) but tools that we have in dimension 3 
(one-to-one mapping with dimer models, determinantal representation) are not available 
to compute them, and therefore the method of [2] cannot be applied. 



4. Proof of Proposition 13.21 

Our proof of Proposition 13. 2| which is the core of this note, uses repeatedly a variant of 
Proposition 3 from [2j that we state below. The result differs from the cited proposition as, 
although we also consider a dynamics in with + boundary condition, we put the extra 
condition that spins inside a sphere of radius S^_i (with / ^ (logr)'^2/2) are constrained 
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to stay — . Nevertheless, the proof in [2] works without any change here. Recall here that 
at the end of Section [2m we fixed C2 = 3/2, ci > 13/2, and cq > ci + 2c2. 



(log'-)'' 



-, consider the (three-dimensional) dynamics in 



Proposition 4.1. Given r and I ^ ^ 
S'^ \ S'^^i with + boundary condition on dS"^ and — on d~ S'^_i (i.e. + on the outside, 
on the inside). Then for any k > 0, uniformly in I ^ (logj-) ^ 

■ 1 



3te 



-r([ogry^ ,r^ 



3x£ Sf\ Sf_^, cj~(t) 



0(r- 



(4.1) 



Note that by monotonicity, the result has to be checked only for / ^ (logr)'^2/2. From 
this proposition, we deduce the following result, that will be of practical use for us. 

Corollary 4.2. Consider the (three-dimensional) dynamics in S^XS"^ on„„ r^co with bound- 



ary condition + on dS^ and 
[L{logLy\8L{logLy'], 



on d 5'^_2(iogL)'=2 ■ 



r-2(logL)=2 

Then for any k > 0, for any r G 



3tG [L(logL)^°,L3] , 3x(^Sl\Sl 



0{L-^) 



V-{logL)'=2 

Proof. Define 

Ti := inf ^ fL(logL)^i+^2 | 3^ ^ ^3 ^ ^3_^^ ^-^^^ 
with the convention that inf = 00. Equation (14. 2p is implied by 

P(rr(iogLr2] ^l3) = o(l-^'), 

because cq > ci + 2c2 . 

We prove by induction on i that for alH G {1, . . . , [(log L)^2"| |^ for all large L, 



(4.2) 

(4.3) 
(4.4) 



For i = 1 one simply uses Proposition 14. II with / = 2(logL)'^2 ^ (log r) 2 



(4.5) 

For the induction 



step, one uses the induction hypothesis as follows 

P(Ti+i ^ L3) ^ P(Ti ^ L^) + P(Ti+i ^ L3 . > L3) 

^ iL-^'^+i) + P(ri+i ^ ; r, > L^). (4.6) 

Now we need to estimate P(t,_|_i ^ ; r-i > Li^). In order to do so we introduce an 
auxiliary dynamics (fT(t))t ^ that coincides with (y~{t) up to time iL([og LY^^'^'^ , and 
which for t ^ ^^(logL)'^^"'"'^^, uses the same clock process and the same coin-flips (see 
Remark 12. ip as {(y{t))t^ i{\ogLYi+''2 but cancels the moves that create a — outside Sl_^. 
Notice that up to time Tj, a~{t) and a{t) coincide and therefore 



P(Ti+i ^ ; r, > L^) 
•(3tG [{i + l). 
^P(3t G [(i + l)L(logL) 



P(3tG [(^ + l)L(logL)'=l+^^L3] ,3xG a,(t) = -; r, > L^) 



?3 ;^ _ ^ ^3^ 

V-(i+i; 

^^^+^^L3],3xG53\5^(,+l), a.(t) = -; 

Vx G S3 \ a-(iL(logL)^2+^i) = +), (4.7) 

where the last line just uses 

{Ti > L^} C {Vx G \ S^_i, a-{iL{log 1)"^+"^) = +}. (4.8) 
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Now let ^ ^ denote the dynamics in S'^_. \ 5'^_2(iog l)=2 with boundary condition 

+ on dSf_^ and — on d~ Sf_2(iogLy2 started from ^ and denote by P the associated 
probabihty. Using the Markov property at time iL{log LY^'^'^^ , one gets that the last line 
of ()4.7p is equal to 

E 



l{Vxe53\e,..-(.L(iogL)=.+^i)=+}P(3i G [L(logL)^^+^^L3 - iL(logL)^^+^^] 
3x G e(m)>r(^^('°^^)"'"HO = 



(4.10) 



^ P(3t G [L{logLr+'^,L^ - iL{logLr+'^] 3x € St^ \ Kit) = -) (4.9) 

where in the first line a~ {iL{log Ly^^'^^) is, with some abuse of notation, considered as 
an element of <?3 \ <?3 . The second line is obtained by monotonicity. 

•='r-i\'='r_2(logL)=2 

Note that, provided that L is large enough, for every values of r and i that we consider, 
((^ - i) - (r - 2(logL)^2)) ^ i(log(r - i)^ 

and 

L(logL)^i+^2 ^ lr(logr)'=i 
lb 

so that one can use Proposition 14. II and get 

P {3t G [L(logL)^l+^^L3 - iL(logL)^i+^V2] 3x G S',_, \ S'^_^^^,^,d- (t) = -) ^ L-^-^ 
and therefore 

P(Ti+i < L3 ; > L3) ^ L-'^-l. (4.11) 
Plugging (|4.1ip into (|4.6p gives (j4.5p for i + 1 and ends the induction step. 

n 

We are now ready to prove Proposition 13.21 The main idea in the proof is to control 
the dynamics by cutting Cl into L three-dimensional slices and to control the evolution 
of each slice by using the results we have for the dynamics in three-dimensions. However 
this stochastic domination is not straightforward and comes from the ad hoc construction 
of the sets C^*\ 

Remark 4.3. The geometric strategy we use to control the evolution of the set of — in Cl 
presents some similarities with the one used by Caputo, Martinelli and Toninelli to prove 
Theorem 4.1 in [1] concerning the mixing time of the dynamics of volume-biased plane 
partitions. 

Recall that we consider the dynamics in the cylinder Cl with boundary condition rjQ 
described in the introduction (+ everywhere except at the top). Define 

Ti := inf|t^iL(log L)^» I 3xGCl\c['\ fj-(t) = -}. (4.12) 

Proposition 13.21 can be expressed as: for all i G {1, . . . , 4L}, for every k and for sufficiently 
large L (how large depending on k) 

P{Ti ^ L^) ^ iL-'^^+^\ (4.13) 
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We prove it by induction on i. 

We start with the case i = 1. In that case, one just has to show that the set of — in the 
first shce of the cyhnder Cl has decreased, i.e. that (recall notation (|3.8p ) 

P{3t G [L{logLr\L%3x G {Sl,o\Sl,i) x {!}, a-{t) = -) < L-^^+^\ (4.14) 

By monotonicity in the boundary condition, the evolution of the spins on Sl^q x {1} for 
the dynamics on Cl with ryg boundary condition dominates the evolution of the spins for 
the dynamics on Sl^ x {1} with + boundary condition on {Sl^ x {0}) U {dSifi x {1}) 
and — on S^^q x {2}. 

One can then check that the latter dynamics corresponds (up to a trivial bijection) 
to the (three-dimensional) dynamics on Sl^ with + boundary condition. Indeed, for 
each X G Sla x {1}, due to boundary condition one has ax+ei{i) = — , (yx-e4,{t) = + 
(where 64 := (0, 0, 0, 1)), so that the effects on the dynamics of the neighbors in the fourth 
direction above and below cancel out. Once one has noticed this, the case i = 1 is an 
immediate consequence of Corollary 14.21 (without the restriction of having spin blocked 
inside the inner-sphere). 

For the induction step, we use a strategy similar to the one used in the proof of Corollary 
14.21 First, notice that 

P(T,+i < L^) ^ P(r, ^ L^) + V{Ti+i ^ r, > 

^ iL-C^+i) +P(ri+i ^L^ ■,T,> L^), (4.15) 

where we used the induction hypothesis in the last inequality. Now we need to estimate 
P(rj_|_i ^ ; Tj > L^). In order to do so we introduce an auxiliary dynamics {&{t))t ^ 0, 
that coincides with {a~(t))t ^ up to time iL{logLy° and which for t ^ iL(logL)'^'' uses 
the same clock process and the same coin-flips (see Remark l2.ip as {(7{t))f ^ jL(iogL)'=o but 

cancels the moves that create a — outside of Note that a{t) and a{t) coincide up to 
time Tj. Therefore 



P(T,+i < l3 ; T, > L^) 

= P{3t G [{i + l)L{logLr,L^] 3x G Cl\4^'\ cT,.(t) = - ; > L^) 
^P(3t G [(i + l)L(log L)^«,L3] 3xGCl\4'+'\ ax{t) = - ; 

\fx G CL\4\(T-{tL{logLr) = +). (4.16) 

Now let (d-~{t))^ ^ g denote the dynamics in c|*^ with — boundary condition on {Sq^i^^^^c2 x 

{L + 1}) n dC^^^ (the top) and + everywhere else. As in the proof of Corollarv 14.21 by the 
Markov property at time iL(logL)'^" and monotonicity the last line of (|4.16p is at most 



3t G [L{\ogLr,L^ - iLilogir] 3x G C^;^ 4*+'^ a-{t) 



(4.17) 



Now notice, that with our particular construction of C^*^ (see Figure [2]) the j-th. slice of 
cj*^ is equal to the j + 1-th slice of C^^'^'^^ for all j ^ f^/^l j so that for every x G cj*^ \ C^^~^'^^ 
one has x it 64 ^ cj*^ \ C^l~^'^\ 
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More precisely one has for every i € {0, . . . ,2L — 2} 



U ((9"5i,(i+2)(iogL)=2 X {1}) U {{dSL,o) X {ri/21 + 1}) . (4.18) 

The same holds for i € {2L — 1,...4L}, but with {dSLfi) x {\i/2] + 1} replaced by 
(552,^(j_2(L-i))) X {L}. Here U denotes the disjoint union. 

Let (cj(t))t ^ denote the dynamics in cj*^ \cj*^^^ with the following boundary condition 
. - on (C« \ 4'+')) + 64 and on (5-5^,(^+2)) x {1}, 

. + on \ 4+')) - 64 and on (dSLfi) x {[^21 + 1} or (5i,(,_2(L^i))) x {L} 

(see Figure [2] for a graphical description of the boundary conditions). Once again, mono- 
tonicity in the boundary condition gives that (j4.17p is at most 

3t e [L(logL)^o,L3] 3x G C« \ 4+'\ ?-(t) = -] . (4.19) 

To finish the proof we just need to prove that the above quantity is less than L^^^^. 

The key point is to notice that the boundary of cj*^ \ C^^^'^'^ is the union of the boundary 
of its slices , i.e. 

L 

d (Cf \ Cf^)) = \Jd [(5,.,(,_20-i)), \ 5i,(.-20--i)+2) J X {j}] . (4.20) 
i=i 

Therefore the evolution of the spins in the different slices are independent. Moreover, 
due to our choice of boundary condition, the evolution of {a{t))t va. the j'-th slice of 

{SL,{i-2{j-l))+ \SL,(i-2{j-l)+2) + ) X {]} 

is the same (up to a trivial bijection) as the evolution of the spins for the dynamics in 

SL,{i-2{j-l))+ \ 5L,(i_2(i-l)+2)+ 

with + boundary conditions on 95/^^(j_2(j-i))+ aiid — on c^~5'/,^(j_2{j-i)+2)+ • 

The reason for that is that for every x G {Si^(j^^2(j-i))+ \SL,{i-2{j-i)+2)+) x {j} the 
boundary condition we have chosen imposes that (Tx-ca = + ™d cj^+eA = — and therefore, 
influence of neighbors along the fourth direction cancels out. For the rest of the boundary, 

imposes + boundary conditions on (c?5/,,{j-2(j-i))+)x{j} and — on (c?~S'j;^^(j_2{j_i)+2)+)x 
{j}. This appears well on figure [2l 

Hence, using again Corollarv 14.21 one gets that for L large enough, for every i S [1,4L], 
j G [1, -Z^], one has 

P[3t G [L(logL)^«,L3] 3x G {SL,(i-2ij-i)u\SL,ii-2U-i)+iu) x {j}, 

a-{t) = -] ^ L-('=+2). (4.21) 
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+ , ,+ 

+++ I I - TJZ, .JIZ- ~ , I +++ 

+++ , I 1 , +++ 

+++ 1 _~ ~ ~ , +++ 

: . zir. . j ++" 

+++ +++ 

Figure 2. Schematic two dimensional view of the boundary condition considered for 
the dynamics > q (for i ^ 2L — 2, i even). Here the full line denotes the border 

of C^'^^', the thick dotted line denotes the border of C^*' and the shadowed region is 
C^'' \C^*^^''. The boundary conditions are also represented: for every site in C^'' \C^'^^\ 
the neighbors along the fourth direction are outside of C['' \C['+^' and the two neighbors 
along that direction are given opposite spins by the boundary conditions. 

Using the union bound and summing over j G {1, . . . , L} one gets that 



^ P 



3t G [L{logLr,L^] 3x G a' (t) = - ^ L-^^+^\ (4.22) 



which finishes the proof. □ 

5. From 4 dimensions to d dimensions 

We briefly indicate how the proof should be modified to work also in higher dimensions. 
One replaces by 

CL,d ■= S2dL(\ogLY2 X {1, . . . ,lY ^- (5.1) 
One considers the dynamics in Ci^d with boundary condition 

• + on (^^2dL(iogL)=2) X {1,...,^}'^"^ and on Uti 52dL(iogL)-2 X Fi where F, := 
{l,...,Lr-ix{0}x{l,...,L}'^-3-i 

• - on U ti 52dL{iogL)^2 X F[ where i^' := {1, . . . , Lf-^ x {L + 1} x {1, . . . , LY-^-\ 
Define: 



'^L,d ■~ '-'2dLlogL'=2-{i-k+2{d-3))+{logL)'=2^ 



C'l!d-=^L,d'~^{^={^1^^2,Z3,Z{d-3)) ^^'^ \ {Z1,Z2,Z3) ^ S'l'^ '''^| 

where z^'^~^^ = (24, . . . , Zd) G Z'^"'^ and | • • • |i denotes the li norm on Z'^^^ 



(5.2) 



\z^'-%:=Y.\^,\. (5.3) 

j=4 

Then one can prove the following generalization of Proposition 13.21 
Proposition 5.1. For every k > 0, for every i G {1, . . . , 2{d — 2)L} one has 

P \3t G [iLilogiy^L^] 3x G Ci,^ \ (0 = "l = 0(1-^). (5.4) 
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The proof is performed in the same way, using induction on i. Note that the proof is 
still based on the three-dimensional result Corollary 14.21 and not on an higher dimensional 
version of it (that we would be unable to prove). Coupling, stochastic comparison and 
identification with the three-dimensional dynamics are valid in that case also. The fact 
that on any intermediate result, the value of k (when one shows a quantity is less than 
L~^) can be chosen arbitrarily large, allows to perform a union bound in any dimension 
without any harm. In the induction step, using coupling and stochastic comparison, one 
is left to study the dynamics in C^l\ \ C^l~d^^ with specific boundary condition. 

Note that in analogy with (j4.18p the boundary of i9(cj*'*^ \ cj*^^^) can be decomposed 
in four parts as follows 



{i+2) 



u 



u 



LA \ '^LA 



X {iC'i-s)}} u/c, (5.5) 



where l^'^ = (1, . . . , 1) G Z*^ ^ and /C denotes what remains of the boundary when the 
three other parts have been taken away. 

Our specific boundary condition is 

Cj) and /C, 



+ onU?=4(Ca\c(:^^ 



• - on Uj=4(CL,d \ ^Id - '='3 J '^-^ \^ ^L,, 

One can check that with these boundary conditions, the different slices 

of the system evolve independently for each z^'^~^^ £ {1, . . . , L}'^^^ and have the same 
evolution as the corresponding three-dimensional dynamics with appropriate boundary 
condition (for j G {4, . . . , d} for every x G ^ =t belong to the boundary and 

our boundary condition ensures that the influence of spins of x it ej cancels out). One can 
apply Corollarv 14.21 and perform a union bound on z^'^~^\ Details are omitted. 
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